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We investigated travelling waves appearing as the primary pattern-forming in¬ 
stability in the nematic Phase 5 (Merck) in the planar geometry in order to 
test the recently developed weak electrolyte model of ac-driven electroconvec¬ 
tion [M. Treiber and L. Kramer, Mol. Cryst. Liq. Cryst 261 , 311 (1995)]. 
Travelling waves are observed over the full conductive range of applied fre¬ 
quencies for three cells of different layer thickness d. We also measured the 
elastic constants, the electric conductivity, and the dielectric constant. The 
other parameters of Phase 5 are known, apart from the (relatively unimpor¬ 
tant) viscosity aq and the two parameters of the weak electrolyte model that 
are proportional to the geometric mean of the mobilities, and the recombina¬ 
tion rate, respectively. Assuming a sufficiently small recombination rate, the 
predicted dependence of the frequency of the travelling waves at onset (Hopf 
frequency) on d and on the applied frequency agreed quantitatively with the 
experiments, essentially without fit parameters. The absolute value of the 
Hopf frequency implies that the geometric mean of the mobilities amounts to 
1.1 x 10 -10 m 2 /(Vs). 

PACS 47.20-k, 47.65+a 
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I. INTRODUCTION 


Electrohydrodynamic convection (EHC) in nematic liquid crystals (NLC) is one of the 
most prominent phenomena employed for the study of pattern formation in anisotropic 
systems []T]-[|- The system consists of a NLC with negative or only mildly positive dielectric 
anisotropy sandwiched between two glass plates with transparent electrodes. When applying 
an ac voltage with the circular frequency ujq and increasing the rms V above a certain 
threshold V c (cu 0 ), the non-convecting basic state becomes unstable to stationary or travelling 
rolls. In this work, we consider the usual planar-homogeneous geometry where the director 
of the NLC is aligned parallel to the plate electrodes. 

The traditionally used hydrodynamic model of EHC originating from Helfrich Q in¬ 
volves equations for the velocity held, the director, and the electric potential (or charge 
density). The NLC is considered as an anisotropic ohmic conductor. The origin of the 
conductivity is not specified. This standard hydrodynamic description M, including the 
three-dimensional formulation for EHC p,|3]U, is referred to as the Standard Model (SM). 
For a general introduction, see, e.g., the book of Blinov |L0|, or the reviews |I[. 

Many features in the conductive range of low ac frequencies are quantitatively described 
by the SM [[|. Examples are the threshold voltage as a function of the ac frequency u; 0 , 
the wave vector of the pattern, which describes the spacing and the angle of the rolls and 
includes the possibility of oblique rolls (where the roll axis is not perpendicular to the 
direction of planar alignment), [||,[| and, in the nonlinear regime, the amplitude of the roll 
pattern and secondary bifurcations leading in particular to (weak) defect turbulence |TT|-[l3|. 
Other features remain unexplained even on a qualitative level. Most notable are travelling 


waves (TW) which have been observed as early as 1978 [Tj|. Later on, they were found by 
various groups in a broad parameter range in the NLCs MBBA ||T5| , [T6|1 , 152 |]T7| , |l8| , and in 
Phase 5 [2|.[T0l| . TWs appear to be generic for relatively thin and clean cells and seem to be 
often favoured by higher external frequencies. 

The recently developed weak electrolyte model (WEM) [^0| predicts a Hopf bifurcation 
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for these conditions, provided the recombination rate of mobile ions (see below) is sufficiently 
low. The WEM equations for the director and velocity fields are those of the SM, but the 
assumption of ohmic conduction is dropped and the conduction properties are explicitly 
modelled by two species of oppositely charged ionic charge carriers. The dynamics of the 
number densities n + (r,t) and n~(r,t ) of the charge carriers is governed by drift relative to 
the fluid under an electric field and by a conventional dissociation-recombination reaction 
between neutral molecules and the ions. The WEM expresses the total space-charge density, 
which already appears in the SM, as the difference of the number densities of the two ionic 
species. In addition, the drift of the charges excites a new field, which can be expressed in 
terms of a variable local conductivity cr±(r,t) perpendicular to the director (’’charge-carrier 
mode”). 

The SM has two relevant time scales: the director-relaxation time 77 = 7 1 d 2 /(A' 1 i 7 r 2 ) 
and the charge-relaxation time r q = eoej_/<T_L. Here e± is the principal value of the dielectric- 
constant tensor perpendicular to the director, 71 is a rotational viscosity, and K n is the 
orientational elasticity for splay distortions. Typically, 77 is 0(1 s), and r q is 0( 10~ 2 s). 
The WEM has two additional time scales ma a recombination time r rec for the relaxation 
towards the equilibrium of the dissociation-recombination reaction where crj_(r, t ) = cyjj 1 , and 
the transition time r* = d 2 /(V/j,) for a charge with the mobility fi = + /rj to traverse the 

cell under the critical applied voltage. The mobilities and /ij are the principal values of 
the mobility tensors perpendicular to the director, see Eq. (||) below. The magnitude of T t 
will turn out to be of the order of 0.1 s whereas r rec > 5 s. 

The SM is recovered in the limits 77/77 >> 1, r t uJo » 1, and 777 \J Ai /» 1- The 
ratio 77/77 is a measure for the relative change of the charge-carrier densities induced by 
the space charges of the Carr-Helfrich mechanism. In our case the third inequality will be 
violated. 

Assuming a recombination rate 77 “/ that is small compared to the circular Hopf fre¬ 
quency uh = 27t f Hl the WEM predicts the Hopf frequency fn to be proportional to 
(^;) 1 / 2 (^l q )- 3 , and depending otherwise only on uq and on the material param- 
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eters of the SM. In experiments on the nematic 152 [|17|, the WEM explains succesfully the 
dependence of the Hopf frequency of the observed TWs on cq’ and on u 0 . Measurements 
on a thicker cell seem also to confirm the predicted d dependence, but no systematic mea¬ 
surements were made. Also, not all of the twelve material parameters of the SM are known 
for 152. Six parameters had to be fitted to the curves of the threshold voltage and of the 
roll angle as a function of cuo- These curves can be calculated, in a good approximation, by 
the SM. 


In this work we present, for the nematic Phase 5 (Merck), tests of the WEM essentially 
without fits of the SM parameters, and we investigate systematically the d dependence of 
the Hopf frequency. 

The threshold V c , the wavevector q c , and the Hopf frequency u>h of the travelling rolls are 
measured, as a function of uq, for three cells of different layer thickness d. We also measured 
the orientational-elastic moduli for splay Ku, twist K 22 , and bend /I33, the equilibrium 
conductivity cy[ q , and the dielectric permittivity e± perpendicular to the director as well as 
their anisotropies ((7® q /cr^ q and e a ). Then, together with known data for the viscosities pi 


all SM parameters of Phase 5 are known apart from the small Leslie coefficient a±. ol\ and 
some remaining uncertainties in the equilibrium conductivity are fixed by fitting V c and q c 
within the SM. Under the assumption of a small recombination rate (as in 152), this enables 
a quantitative test of the dependence of uh on cu 0 . The three layer thicknesses allow a test 
of the predicted d dependence which is independent of the material parameters. Only the 
overall magnitude of ujh is fitted once with (p^pj)) 1 / 2 entering as a multiplicative factor into 
the WEM expression. The predicted and the measured Hopf frequency agree, within the 
experimental errors, for all applied frequencies and all three cells, while the magnitude of 
u>h varies by a factor of more than 50 (see Fig. 4). 


I 11 Section II, we describe the measurements, and in Section III, the relevant predictions 
of the WEM. Theory and experiment are compared in Section IV. Section V concludes with 
a discussion and points at some observed nonlinear effects which may be contained in the 
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WEM. 


II. THE EXPERIMENT 


We used three cells labelled C, D, and E, respectively with areas of the rectangular plate 
electrodes between 10 x 14mm 2 and 14 x 16mm 2 defining the directions x and y. Planar 
alignment along x was obtained with a rubbed polyimide him which coated the transparent 
electrodes. The cell thickness in the z direction was measured on the empty cells by an 
infrared spectrophotometer. Due to the large spotsize of the illuminating light beam this 
method yields an average thickness of dc = 24.3 pm, do = 14.7pm, and d® = 29.3 pm. 
Interference fringes indicated that the local thickness varied by ±2 pm around these values 
throughout the cell due to deformation of the plates at sealing. 

The capacitance and conductance of the cells were measured by a capacitance bridge 
shortly after filling with the nematic mixture Phase 5 (Merck) as well as after finishing the 
experiments (6 weeks later). Comparing these data with the capacitance of the empty cells 
and neglecting stray capacitances, the dielectric permittivity e± and the electric conductivity 
cr^_ q at equilibrium could be obtained. These data are summarized in Table I. 

The cells prepared for the measurements are too thin to determine the exact values of 
the anisotropies of the dielectric permittivity and the electric conductivity. Therefore we 
measured these parameters in a magnetic field using a thick (> 400 pm) cell calibrated with 
benzene. Table II shows the resulting data as a function of temperature. 

The experimental apparatus for studying travelling waves consisted of a temperature- 
control stage with a stability of 0.01 K, electronics for applying the ac voltage, and a video 
camera with a computer-controlled imaging system. All measurements were made at (30 ± 
0.01)°C. 

The convection rolls were visualized by the shadowgraph technique |22|]. The cell was 
illuminated by polarized light with polarization along the director, and the resulting shad¬ 
owgraph signal was monitored by a video camera attached to a microscope using a lOx 
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objective. The images correspond to an area of 257pm x 271pm of the cell. The microscope 
was focused to the top of the cell to allow for a dominant linear contribution of the signal 
22 |, where the wavelength of the shadowgraph signal coincides with that of the director 


distortion. Second harmonics due to the nonlinearities of the imaging technique could be 
seen well above threshold. The video camera was positioned to get an image of vertical rolls 
in the normal-roll regime where q y = 0, i.e., the rows of the pixels were parallel to the wave 
vector of the pattern in this regime. A few measurements were carried out in the regime 
of oblique rolls, so the ac frequency, where the transition between normal and oblique rolls 
occurs (Lifshitz point), could be estimated. 

In the experimental runs, shadowgraph images were recorded for each cell at various 
ac frequencies in the conductive regime. The measuring program for one frequency was as 
follows. At first, the threshold voltage V c was roughly located by an automatic iterative 
search. Then, a background image was recorded much below threshold. This background 
was subtracted on-line from all later images to remove intensity inhomogeneities due to 
illumination and surface irregularities. To obtain an improved value of the threshold voltage, 

the voltages were scanned in steps of Ae « 0.002 from a value above threshold (e ~ 0.02) 

_ 2 _ 2 

to a value below (e ~ —0.02), and vice versa. Here, the control parameter e = V /V c — 1 
denotes the deviation from the threshold. Image recording followed the change of the applied 
voltage with a delay of about 60-90 seconds. 

The wavenumber in the x direction q x and the frequency of the TWs at a given value of 
e were determined by Fourier transforms of a space-time image. To obtain the space-time 
image, a particular row was selected in the central portion of a series of 256 images taken 
at consecutive time steps and this row was copied from all these images into consecutive 
rows of a new image yielding a 256 by 256 space-time plot. Figure 1 shows such a plot for 
three values of e. Figure la indicates that small-amplitude fluctuating patterns could be 
observed slightly below threshold. Travelling waves appear in these images as stripes tilted 
toward the direction of travelling. Both left and right travelling directions can be located. 
The horizontal periodicity in these space-time images equals to 2n/q x , while the vertical one 
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corresponds to the period 1// of the travelling waves. 

To increase the resolution, the Fourier transformation was carried out on 2048 points 
with the first 256 points taken from the image (after subtracting the average intensity), 
and the remaining points padded with zeroes. For each row, the power spectrum P(q x ) of 
the spatial frequencies (i.e., the magnitude of the squared Fourier amplitude) was calcu¬ 
lated and averaged over all rows. The wavenumber q x of the pattern was determined by a 
weighted average over all points in the neighbourhood of the peaks of the power spectrum 
corresponding to the first and second spatial harmonics, 

_ _ 'P'q x P(qx) + \T, n q x P{q x ) 

Qx ~ E ’P(q x ) + E ”P(q x ) 

Here, P(q x ) is the average of the power spectrum over all rows attributed to the wavenumber 
q x , S' and E” sums over all points with a maximum distance of eight points from the peaks of 
the first and second harmonics, respectively. The Hopf frequency is determined by analogous 
transformations and averaging procedures for the columns. 

The denominator of Eq. (,!]), 

X'P(q x ) + X”P(q x ) :=A 2 , (2) 


can be regarded as a measure of the averaged modulation intensity of the shadowgraph 
signal. 

For small distortions, the modulation amplitude A is proportional to the amplitude of 


the director distortions |22|. The amplitude of the director distortions is expected to scale 
with e 1 / 2 for continuous bifurcations, at least, if subcritical fluctuations play no role. For 

e large enough to neglect subcritical fluctuations, but small enough to remain within the 

—2 

linear regime of the shadowgraph technique, one has A oc e. Figure 2 shows a typical 

_2 

example of A vs. e for cell D at f 0 = 90 Hz corresponding to u 0 r g = 2.06. The threshold 

_2 

was defined by extrapolating the first roughly linear rise of A to zero (dashed line). For 
e > 0.09 one enters the nonlinear regime of the shadowgraph, which is indicated by the 
appearance of a second harmonic peak in the Fourier spectrum. The scatter in the data of 
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Fig. 2 for e > 0.011 is a result of the disorder appearing in the travelling pattern. There are 
structureless (nonperiodic) areas separating the left and right travelling waves with random 
space and time distribution which reduce the averaged intensity modulation accidentally. 

The wavenumbers and the frequencies were calculated according to Eq.(l) as a function 
of e. The values q cx and u>h at threshold (e = 0) are obtained by fitting a straight line to 
points in the range —0.01 < e < 0.02. To monitor the change of the conductivity during 
the measurements, the transition frequency u > c d to dielectric rolls (see Sec. IV below) was 
determined at the start and at the end of the experimental run for each cell. u; c d at the time 
of the measurement was determined from these values by a linear interpolation. 

Despite careful preparation of the cells, there was a spatial variation of the thresholds, i.e., 
on increasing V from subcritical values, convection did not fill the whole cell simultaneously. 
To obtain the most homogeneous image possible, a local minimum or maximum of the 
threshold was selected in the cell for the monitored area of the measurements. In cell 
C, a local minimum was chosen (convection sets in first), while in the cells D and E, a 
local maximum was chosen (convection sets in last; rolls appeared first at the edges of the 
electrode). 


Elastic moduli of the substance were determined by measuring the Freedericksz thresh¬ 
olds in various geometries [|7y23|]. The magnetic field induced splay and twist Freedericksz 
transitions could be observed on the same planar cells used for studying the travelling waves. 
For magnetic and electric field induced bend transitions homeotropic cells were prepared. 
The field induced deformation was detected by optical methods (birefringence for splay and 
bend, depolarization for twist transitions). As the anisotropy of the magnetic susceptibility 
is unknown for Phase 5 the measurements of the threshold magnetic field yielded only a 
ratio of the elastic moduli, = 0.47 ; Aaa _ i 29. Measuring the threshold voltage in 
homeotropic cells we obtained A '33 = 12.7 x lO” 12 ^. The error of these measurements, 


which are summarized in Table III, is ±10%. 





TABLES 

TABLE I. Conductivity, dielectric permittivity, and related quantities in the cells used for the 
EHC experiments (T = 30°C). 


Quantity 

Cell C 

Cell D 

Cell E 

Comments 

Average thickness d [pm i] 

24.3 

14.7 

29.3 

There is a ±2/xm deviation 
from d throughout the cell 

Local thickness d [/zm] 

23.5 

13.4 

27.5 

Deduced from the wavevector (see Sec. IV) 

Conductivity (7^ q 

8.5 

16.4 

7.3 

in fresh cell 

[10 -9 (fhn) _1 ] 

8.6 

14.7 

11.1 

6 weeks later 

Permittivity e± 

5.34 

5.38 

5.24 

in fresh cell 


5.30 

5.39 

5.19 

6 weeks later 

Range of the transition 
frequency w c d [Hz] 

84-90 

138 - 177 

112 - 130 

minimum and maximum values 
during measurements 

<7 [10 _9 (Dm) _1 ] 

7.0 

13.0 

8.0 

Average conductivity obtained from 
the measured td c d (see Sec. IV) 


TABLE II. Dielectric permittivities and electric conductivities of Phase 5 as a function of 
temperature for a thick cell (400/mr) at 1 kHz 


T(° C) e_i_ 


£7^(10 9 (ffrn) x ) o\j q (10 9 (flm) x ) <7® q (10 9 (fim) x ) cr^/cr” 1 


20 

5.336 

5.100 

-0.236 

2.89 

4.98 

2.09 

0.723 

22 

5.319 

5.086 

-0.233 

3.07 

5.35 

2.28 

0.743 

23 

5.305 

5.084 

-0.221 

3.29 

5.65 

2.36 

0.717 

25 

5.286 

5.079 

-0.207 

3.66 

6.44 

2.78 

0.760 

27 

5.256 

5.058 

-0.198 

4.25 

7.38 

3.13 

0.736 

30 

5.217 

5.033 

-0.184 

5.16 

8.73 

3.57 

0.692 

32 

5.185 

5.013 

-0.172 

5.89 

10.17 

4.28 

0.727 

35 

5.159 

4.994 

-0.165 

6.69 

11.22 

4.53 

0.677 
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FIGURES 



£=-0.001 £ = 0.005 £=0.017 


FIG. 1. Space-time images of the travelling waves in cell D at / = 90 Hz. The contrast was 
enhanced by a factor of 5.3 for all three images, (a) Below threshold (e = —0.001). Subcritical 
fluctuations can be seen, (b) Sligthly above threshold (e = 0.005). (c) Further above threshold 
(e = 0.017). Travelling waves in both directions are observable. 



applied voltage for cell D at / = 90 Hz. The threshold V c is defined by extrapolating the first 
roughly linear rise of I to zero (dashed line). The three squares denote the points corresponding 
to the space-time images in Fig. 1. 
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III. PREDICTIONS OF THE WEAK ELECTROLYTE MODEL 


In the WEM, the conductivity of the NLC is described by the drift of two species of oppo¬ 
sitely charged freely mobile ions originating from a dissociation-recombination reaction from 
neutral molecules. The ionizable molecules are not specified. They can be (uncontrolled) 
impurities (which is assumed in most experiments on MBBA and Phase 5 and presumably 
accounts for the strong variations of the conductivity), or ionizable dopants (which is more 
controlled and realized, e.g., in 152 [|l7l , |I§[l ). 

The WEM assumes that only a small fraction of the impurities is dissociated into ionic 
charge carriers (weak electrolyte). This is not an essential restriction but it allows to reduce 
the number of material parameters not contained in the SM. For typical values for the mo¬ 
bilities [10~ 10 m 2 /(Vs)] and conductivities [10” 8 (flm) _1 ], the number density of both species 
of ions is about 3 x 10 20 m -3 . For a molecular weight of the impurities of 100, this corre¬ 
sponds to a concentration of 0.05 ppm. One has a weak electrolyte, if the number density of 
the impurity molecules is sufficiently higher. In the experiments employing 152 mentioned 
above, this is certainly fulfilled [^4|. In the present work, the weak electrolyte assumption 
is consistent with the roughly exponential increase of the conductivities with temperature 
(Table II) corresponding to a dissociation energy of about 0.4 eV. 

The drift velocities of the oppositely charged ions (superscript + and —) relative to the 
fluid under the influence of the electric field E are given by 

drift = [1% = (if (Sij + (T^niUj) , (3) 

defining the uniaxially anisotropic mobility tensors for the two ions. Here, Sy is the Kro- 
necker symbol, rii and rij are components of the director, and <r' a = cr^/cr” 1 is the relative 
anisotropy of the conductivity. The local conductivity perpendicular to the director (which 
is constant in the SM), is given in terms of the number densities n ± . For charges ±e, one 
obtains 

ct_l (r, t ) = e (ji\n + + pi . (4) 
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The basic equations of the WEM are given in Ref. [^0| . They are obtained by transforming 
the rate equations for the ionic number densities n + and n - into an equation for cr±(r,t) 
describing the dynamics of the drift and of the dissociation-recombination reaction, and 
a conservation equation for the charge density p = e{n + + n ~) together with the Poisson 
equation. The balance equations for the director components and the generalized Navier- 
Stokes equations for the fluid velocities are those of the SM. 

The non-convecting basic state is nearly homogeneous, which simplifies the linear stabil¬ 
ity calculations. A one-mode Galerkin approximation for the linearized fields yields quanti¬ 


tative results which typically deviate by less than 3% from the fully numerical solution [24 


In this approximation, the dependence of each field on z and t is represented by one test 
function satisfying the boundary conditions and the symmetries under reflection and time 


translation applying in the conductive regime |20|. This one-mode approximation enables 
analytic formulas for the Hopf frequency, and for the threshold shift Ae = ( V c /V ^ M ) 2 — 1 
of the WEM with respect to the SM. For details in the normal-roll regime, see Ref. |[20|| ; 
the generalization to oblique rolls can be found in Ref. [T7[. The threshold shift is always 
smaller than 1% for the experiments of this paper, so the thresholds, and also the critical 
wavevectors, are calculated using the SM. 

For CiiTd « VWoej_/(-Rii 7 r 2 ) which is fulfilled for all measurements, the analytic expres¬ 
sion for to H = 277 ffj is 


u>h — AWl — 


[Tree ^) 2 ’ 


U = 77 C 


V:e 0 e ± 


d 3 ( 1 + u/ 2 ) \ 7 iU^ q 


M±M± 


(5) 

( 6 ) 


Here, C' = \/o r a C is a dimensionless factor of order unity with C given for normal rolls in Ref. 
HI ( see also H), and a/ = LO 0 r q (l+qp+e' a q? x )/{l+ <?' 2 +< q? x ) with e' a = e a /e ± , q = (q x , q y ), 
and q' = qd/ir. In the experiments described here, one has always 0.45 < C < 0.75 and 
0.60 0J()T q < CO 1 < 0.68 (jJftTq. 

The main ingredients of Eqs. (|5|) and (^|) can be seen by considering the two feed- 
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back mechanisms resulting from the WEM. The primary instability mechanism is the Carr- 
Hclfrich mechanism of the SM. Under an applied held, an initial director bend leads via 
the anisotropy of the conductivity to a charge accumulation and to an electric volume force 
driving the fluid motion. The resulting velocity gradients couple back to the director via 
the rotational viscosities. The time scale of this mechanism is ~ t ( j. The gradients of the 
charge accumulation, however, also excite the charge-carrier mode whose feedback tends to 
decrease the charge. The time scale of this second stabilizing drift mechanism is ~ T t I T q- 
As in binary mixtures (see, e.g., ||2(|), the interplay of the two mechanisms can lead to os¬ 
cillations with a (Hopf) frequency proportional to the geometric mean of the two processes 
involved. This gives u ~ , which is Eq. (^]), up to a dimensionless prefactor. 

The drift mechanism is only effective if the charge-carrier mode can build up sufficiently, 
i.e., for a small recombination rate. This leads to the square root in Eq. (^) and to the 
condition t~ 1 < u) for the occurrence of a Hopf bifurcation to TWs. The condition u> = r~l 
gives the codimension-two point separating the regime of stationary rolls (r”* > uj) from 
that of TWs. This condition scales with d~ 3 , (cdj) 1 ) -1 / 2 , and V 2 c /{ 1 + u/ 2 ), hence a Hopf 
bifurcation is favoured by thin cells, clean cells (low oy_ q ), and by high values of the applied 
frequencies in case of materials with e a < 0 (for e a < 0, the factor V c /(1 + u/ 2 ) increases 
with Uo and contains the main part of the dependence on cno). 

In contrast to former experiments on Phase 5 j2]JUj, we observed TWs for all applied 
frequencies in the three cells. This is presumably due to a smaller concentration of impurities. 
The Hopf condition is most restrictive for the thickest cell at the lowest frequencies. The 
experimental points for this cell (squares in Fig. 3a) suggest that one is still far away from 
the codimension-two point. With Eq. (^), this means that a >h ~ u), i.e., l/r rec , can be set 
to zero in the calculations. 

In this case, the WEM predictions contain (/i^/i])) 1 / 2 as the only non-SM parameter, 
which appears in ujh ~ d), as a simple multiplicative factor. The WEM predicts that u) is 
proportional to d~ 3 (a[[ q ) -1 / 2 and increases with the ac frequency. In the following section, 
these predictions will be compared with the experiment. 
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IV. COMPARISON BETWEEN EXPERIMENT AND THEORY 


The WEM gives V c , q c , and the Hopf frequency c vh, as a function of the applied frequency 
ujq r q normalized to the inverse of the charge relaxation time T q . The values for V c and q c differ 
by less than 1% from those of the SM. In a good approximation, V c {ojqt^) and q c (u>oT q )d/7T 
do not depend on d or on cr” 1 [j3|], the only parameters that are significantly different in 
the three cells. All experimental points for V c and q c should lie on one curve, if one plots 
V c and q c d/7T against a; 0 T g . Similarly, the WEM requires that the experimental points for 
cuHd 3 (a e ^) 1 ^ 2 should lie on one line. To verify this, requires knowledge of r q oc 1 /a c ^ for each 
experimental point. 

The conductivity turned out to vary somewhat irregularly during the experimental runs 
(Table I). The relative change of the conductivities during the measurements can be assessed 
by determining, for each measurement of the set V c , q c , and ujh, the transition frequency u; c d 
where the rolls of the conductive regime cross over to dielectric rolls ||T[, which can be easily 
identified. This frequency is related to the ’’cutoff frequency” u; cuto ff where the threshold of 
the conductive regime diverges in the one-mode approximation. The ratio a; c d/ct; cu toff depends 
on d, but is constant for one cell (u; cuto ff/u; c d = 0.68- • -0.74 for the three cells). Since the 
cutoff frequency in the SM is proportional to a e _^ (^cutoff = 5.57T' 1 for the parameters of 
Table III), the actual conductivity can be written as 


eq _ 
± ~ 


^cd 

^cd 


O’, 


(7) 


where uJ c d is defined as the average of the transition frequencies over all measurements in 
one cell (C,D, or E) and the ’’average conductivity” a must be determined separately. 

Figure 3a shows the resulting threshold voltages (from the SM) as a function of u 0 r q = 
‘AAoC-l/'W' 1 with <7 l ( q from Eq. (0) (continuous curve) together with the experimental results. 
For each of the cells, a was fitted to the threshold curves of the SM with the other parameters 
taken from Table III in Appendix A. Of course, the resulting value of a must be consistent 
with the independent conductivity measurements (Table I) taken some time before and after 
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the measurements of the patterns. Table I shows that a deviates by less than 25 % from 
these measurements, which is within the variations of the conductivities as estimated from 
the variations of cu c d- The unknown viscosity aq was adjusted as well. In the normal-roll 
regime (all but the measurements corresponding to the three leftmost triangles are in this 
regime), the influence of aq on the threshold and on u>h is negligible. In contrast, the 
tendency to oblique rolls is increased by increasing aq; so aq was fitted to give the correct 
frequency for the transition from normal to oblique rolls, the Lifshitz point (in the SM). 

Figure 3b shows the normalized x-component q^d/'K of the wavevector. Taking for d the 
average thicknesses measured on the empty cells (d, Sec. II) gives a systematic deviation 
(5% - 10 %) of the experimental points from the curves from the SM. One can attribute this 
to the fact that the actual thickness at the position of the measurement may differ from d 
(see Sec. II). The systematic deviations can mostly be removed by using slightly reduced 
thickness values dc = 23.5 /im, do = 13.4/un, and d# = 27.5 /nn(Table I). This correction 
(a reduction by 0.8 /mi for cell C, by 1.4 /mi for cell D, and 1.8 /mi for cell E) is within 
the limits of deviation indicated by the interference fringes (±2 /mi). Hence for the plots 
presented in this paper we assumed these reduced thickness values. The agreement for these 
thicknesses is reasonable, although the theoretical frequency dependence (from the SM) is 
somewhat too steep at larger uor q . The kink in the theoretical curve at uoT q = 0.8 is due to 
the Lifshitz point. 

In Figure 4, we show the main result of this paper, namely the measured Hopf frequencies 
together with the WEM predictions, Eq. (^]), under the assumption of a long recombination 
time where loh = Co. The overall magnitude of loh fixes the mobilities to 

\J /tjj/tJ = 1.1 x 10 _lo m 2 /(Vs), (8) 

but nothing else can be fitted. Since in the normal-roll regime aq plays a minor role for the 
dependence on u>o, the d dependence and the dependence on lo 0 provided by the theory are 
in effect free from unknown or adjustable parameters. 

The measurements for the thickest cell (squares) and the thinnest cell (triangles) agree 
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with the WEM curve nearly everywhere within the experimental error of about 25%; the 
Hopf frequency of cell C (circles) is systematically higher (« 25%) than the predicted one. 
This has to be compared with the wide range covered by the actual Hopf frequencies. In 
the experiments, the Hopf frequency varies from about 0.4s -1 (cell E at uJoT g = 0.8) to more 
than 22s -1 (cell D at u)oT q = 4.0), i.e., by a factor of more than 50. For a given external 
frequency, the Hopf frequency of cell D is about seven times that of cell E. 
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Normalized frequency co 0 x q 

FIG. 3. (a) Threshold voltage, and (b) wavenumber q cx as a function of uoT q . The conductivity 
cr^ 1 in T q is calculated individually for each point from a; c d, see the main text. The experimental 
points of the three cells (circles, squares, and triangles) should lie on one curve for V c and q cx d/ir , 
respectively. The curves are the SM predictions for the material parameters of Phase 5 given in 
Appendix A with the viscosity a\ fitted to —O. 3571 . The WEM predictions for V c and q cx cannot 
be distinguished from that of the SM on this scale. 
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FIG. 4. The Hopf frequency of the travelling rolls multiplied by (cj^) 1 / 2 ^ 3 , is shown as a 
function of the normalized applied frequency ojqTq. The circles, squares, and triangles denote the 
experimental results for the three cells. The curve is the universal WEM prediction for the material 
parameters of Phase 5 (Appendix A) and = 1.1 x 10~ 10 m 2 /(Vs). For ujQT q = 0, this curve 

corresponds to ojh = 0.67 s _1 , 2.65 s _1 , and 0.37 s -1 , for the cells C, D, and E, respectively. 


V. CONCLUSION 

Our results represent the most stringent test for the WEM model so far. The frequency 
dependence of the threshold V c and, to less extent, the critical wavevector qr c , involving in 
essence only SM quantities, has been tested previously fairly well in the conductive range 
Also Qii was fitted to obtain the correct Lifshitz point, but this is of little influence 
in the normal-roll range considered here mainly. The Hopf frequency c Oh predicted by the 
WEM deviates from the measured ones by at most ~ 25% in a range of parameters where 
u>h varies by a factor of about 50. This was achieved by fitting only the parameter I 

appearing as an overall factor. 
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The value of \J given in Eq.(|8]) compares reasonably well with a number of measure¬ 
ments of the quantity /i = in different materials, namely (in units of 10 -lo m 2 /(Vs)) 

fi = 0.18 |28|], 1 [fH| , 3.71 m , 1--10 |3l|, 10 p7[| . Little appears to be known about the ratio 
/Ll/mI- 

Our results indicate that the recombination time r rec is well above 1/loh even for the 
smallest Hopf frequency measured, which means T rec is larger than about 5s. Such long 
times may seem surprizing, however, literature values for MBBA range from 10 _3 s |27J to 
2.7 x 10 4 s E3]. Thus in our experiments one is well away from the codimension-2 point 


that separates the stationary regime from the Hopf regime. Let us point out that a weakly 
nonlinear analysis predicts the bifurcation to be typically subcritical near the codimension-2 


point on the stationary side |24|. This appears to explain the slight hysteresis observed in 


some high-sensitivity experiments ||T7| , |32| , h6| . The fact that a hysteresis is sometimes also 
observed on the Hopf bifurcation side |3^|T6|] is not disturbing because there is evidence that 


this is due to higher-order effects in the amplitude |34|. The theoretical analysis is consistent 
with this fact. 

For further testing of the WEM an independent determination of both r rec and the 
mobilities is desirable. This could be done by measuring the transient current response 
to suitable voltage signals in cells with well defined boundary conditions (e.g., blocking 
electrodes) P3|,|2D|,|28[ . 
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APPENDIX A: MATERIAL PARAMETERS FOR PHASE 5 AT 30 °C 


The parameters K n , K 22 , K 33 , e±, e Q , cr® q , and cr® q were measured in this work. The ac¬ 
tual value of (T^ q varied considerably and was determined individually for each measurement 
by Eq. ('7]) with <7 from Table I. There are published data for the dielectric anisotropy 
which differ somewhat, but these data seem to be rather inprecise. The Leslie viscosities 
ay to olq were determined from the viscosities r/i, 772 , rj 3 , and 71 of Ref. |H|] by the relations 
«2 = (V 2 ~Vi- 7i)/2, «3 = {V 2 + Vi~ 7i)/2, «4 = 2 / 73 , a 6 = (71 + 3 ?/ 2 - 4t ? 3 - 7 i)/ 2 , and 
the Onsager relation a .5 = a§ — a .2 — 01 3 . 
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TABLE III. Material parameters for Phase 5 


Parameter Unit Value at 30 °C Source and comments 

An 10- 12 N 9.8 Sec. II 

X 22 ” 4.6 

A 33 ” 12.7 

(Ja/v^ - 0.7 Table II 


e_L 

ai 

0L2 
a 3 

Oi\ 

«5 


10 3 Ns/rn 2 


11 


5.19-•• 5.39 


-0.184 


-39 


-109.3 

1.5 

56.3 

82.9 


Tables I, II 

Table II; see also Ref. [35| 


Fit to the Lifshitz point 


Kneppe, et al, Ref. [21] 


11 


11 


-24.9 



10 10 m 2 /(Vs) 

1.1 

Fit to ojh 

T~^ 

1 rec 

s" 1 

< 0.2 

Hopf bifurcation in all cells 
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